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Abstract. This papers underscores the intimate connection between the quan- 
tum walks generated by certain spin chain Hamiltonians and classical birth 
and death processes. It is observed that transition amplitudes between single 
excitation states of the spin chains have an expression in terms of orthogonal 
polynomials which is analogous to the Karlin-McGregor representation formula 
of the transition probability functions for classes of birth and death processes. 
As an application, we present a characterization of spin systems for which the 
probability to return to the point of origin at some time is 1 or almost 1. 



1. Introduction 

This paper aims to point out the relation between birth and death processes 
(BDP) and quantum walks generated by spin chain Hamiltonians. BDPs are 
stochastic processes [T7] that are special cases of continuous time random walks 
(CTRW) on a graph. In CTRWs, the steps on a graph only take place between 
vertices related by links and occur at time intervals that are sampled from an expo- 
nential probability distribution. A CTRW is characterized by a generator matrix 
with the jumping rates from one vertex to another as entries. The central object 
of interest arc the transition probabilities as a function of time. 

The continuous time quantum walk (CTQW) associated to a continuous time 
random walk is defined [15], [12], [28] through a quantal Hilbert space represen- 
tation: vertices correspond to state vectors and transitions are governed by the 
quantum mechanics ensuing from taking the (restriction of the) Hamiltonian on 
that set of states to be the (Hermitian form) of the generator matrix of the CTRW. 

In this CTQW framework, the key quantities are the transition amplitudes whose 
unitary time evolutions are thus prescribed by this Hamiltonian. As stressed in 
[12) , the classical transition probabilities and the quantal transition amplitudes are 
basically related by analytic continuation through the familiar t O- it substitution. 

The Markov chains resulting from CTRWs on linear graphs are known as birth 
and death processes [T7]- In these cases, the vertices can be labeled by the non- 
negative integers and transitions only take place between neighboring sites. BDPs 
have numerous applications, for instance in modeling biological populations and 
in queuing theory. Their seminal analysis by Karlin and McGregor [25], [26] rests 
on the fact that the generator matrix is tridiagonal, allowing for the transition 
probabilities to be expressed in terms of associated orthogonal polynomials. (For a 
review see [23"]V 

Karlin and McGregor adapted the ideas that had been used in the case of diffu- 
sion processes by Feller [16] and McKean [32] to study birth and death processes 
on the nonnegative integers in the cases of discrete [27] and continuous [25], [26] 
time parameter. This goes along the spectral theory of self-adjoint operators and 
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orthogonal polynomials play a crucial role. The earlier work featured a second order 
differential operator and now we are dealing with a tridiagonal matrix. 
These ideas have been extended in several directions. 

In [14], [19] one finds two independent developments to deal with the case of 
quasi birth and death processes, featuring block tridiagonal matrices and Krein's 
theory of matrix valued OPs. A special case of this already arises in [27] when the 
authors attempt to deal with the case of the integers, see the last section of that 
paper. 

A different extension is first considered in [7J for the case of a discrete time 
QW either on the integers or the nonnegative integers. In this case the CMV 
representation of the unitary matrix U gives the link to (Laurent) OPs on the unit 
circle. These ideas have been elaborated in [8], [20] . [21] . [22] . 

The present paper can be seen as an effort to use the Karlin-McGregor method 
in the case of quantum walks with continuous time. 

In the context of discrete time one knows how to make a quantum walk out 
of any unitary matrix U given in its CMV form. On the other hand if one is 
given an arbitrary Jacobi matrix it is not always possible to make a BDP out of 
it; extra conditions on the Jacobi matrix are required. We will see later that a 
similar situation arises in the case of continuous time parameter. We will discuss 
an example of a tridiagonal matrix which fails to give a CTRW but gives rise to a 
QW. 

It may be appropriate to note that the similarities as well as the differences 
between quantum walks with discrete or continuous time is an important subject. 
We refer the reader to the recent paper [11] which by extending a quantization 
procedure for discrete time Markov chains due to Szegedy |33j establishes a corre- 
spondence between continuous and discrete time QWs on arbitrary graphs. In the 
classical case the use of spectral methods, i.e. Karlin McGregor's formulation is in 
principle restricted to reversible Markov chains. 

A hrst point of this article is to underscore the fact the CTQWs corresponding 
to BDPs can be viewed as single excitation dynamics of spin chains. These systems, 
among which one has the Heisenberg spin chain, provide fundamental models and 
are at the heart of significant advances in mathematical physics. They have recently 
been introduced in the realm of quantum information in the design of (perfect) 
quantum wires [6] , [13] . 

In particular, it has been shown [3], [24] that with properly engineered couplings, 
spin chains could effect perfect state transfer, i.e. transport a qubit from one loca- 
tion to another with probability one, without the need of external control. For such 
spin chains with couplings that vary from site to site, the hopping of single excita- 
tions amounts to CTQWs that parallel BDPs. The analogy carries mathematically 
[34] , [9] and here again orthogonal polynomials provide a representation of the time 
dependent transition amplitudes which is similar to the Karlin-McGregor formula 
for the BDP transition probabilities. 

A further point of the paper is to identify the conditions for perfect or almost 
perfect return. This respectively refers to cases where there is probability 1 or 
almost 1 to return to the initial state of the quantum walk after some time. 

The remainder of this paper is organized as follows. In the next section, we 
review the essentials of the Karlin-McGregor treatment of BDPs in a way that 
lends itself to the comparison we seek to make. In section 3, we indicate precisely 
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how the CTQWs corresponding to BDPs are connected to the quantum dynamics 
of spin chains. In Section 4, the natural issue of return is addressed. It is studied 
using the spectral representation and characteristic functions (Fourier transform) 
of the weight distribution of the associated OPs. It will be seen that perfect return 
requires the weight function to be purely discrete on a uniform lattice, while a pure 
point spectrum ensures almost perfect transfer. 

Section 5 provides a number of examples of systems exhibiting perfect return. It 
is noted in particular that the uniform Heisenberg XX chain does not satisfy the 
condition for almost perfect return. 



2. The essentials of birth and death processes 

Birth and death processes are stationary Markov processes whose state space 
is the set of non-negative integers. For BDPs, the transition probability functions 
form a matrix P whose elements Ptj (t), i, j = 0, 1, 2, . . . (defined as the probabilities 
that the system evolves from state i to state j in time t) are required to satisfy the 
following conditions for infinitesimal time (t — > + ): 

(2.1a) P iJ+1 (t) = kit + o(t) 

(2.1b) P iti (jt) = l-(Xi+iM)t + o(t) 

(2.1c) Pi,i-i{t) =fj,it + o(t) 

and Pij{t) = o(t) for \i — j\ > 1. Obviously Pij(0) = &y. Transitions thus only 
occur between nearest neighbors. The coefficients A.; and fj,i assumed to satisfy 
A; > 0, /Ltj+i > for i > and > are referred to as the birth and death 
rates at state i, respectively. By hypothesis the conditional probability of the 
process up to time t depends only on the state of the process at time t, this implies 
that P i:j (s + t) = T,k P ik{s)Pk 3 (t), i.e. P{s + t) = P(s)P{t). This fact, together 
with (|2.1a|) - (|2.1cp readily leads to the (forward) Chapman-Kolmogorov differential 
equation which is obeyed by the transition probabilities. Indeed one has 

P tJ (t + 5t) = J2 Pik{t)P kj (5t) = 

k 

(2.2) P,, J -iWF r i J (fl) + P iJ+ i(t)P j+ i d (5t) + Pi^PniSt) + 0(5t) 
from where we get 

(2.3) = Xj-iPij-iit) + H+1 Pi, j+1 {t) - (A, + H )P l3 (t) 
In matrix form, (|2.3I) can be written 

(2-4) f = ^< 

where 

-(Ao+Mo) A 

(2.5) A= I V\ -(Ai +Mi) Ai . 
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is hence the generator of the CTRW. In order to solve (|2.3|) (or (|2.4|) ). following 
Karlin and McGregor [25] . we introduce the set of orthogonal polynomials Q n (x) 
satisfying the 3-term recurrence relation 

(2.6) - xQj(x) = ^jQj^i(x) - (Xj + fij)Qj(x) + XjQj^i(x), 
with Q-i = and Qo = 1. In compact form 

(2.7) - xQ(x) = AQ. 

It is convenient to introduce the set of polynomials {x^^)} related as follows to 
the set {Qi(x)} 



(2.8) 



fa) = J2 Ut ^^ :> 



with 
(2.9) 
and 

(2.10) 

They satisfy 
(2.11) 



Ui 



(-l)V 



1/2 



A Ai . . . A,;_i . 
7T = 1, 7Tj = , I > 0. 

A*iM2 ■ • ■ Mi 



with J the symmetric Jacobi matrix 



(2.12) 

where 
(2.13a) 

(2.13b) 



J = -UAlf- 1 



{B 










..A 


Ji 


B x 


■h 










B 2 


J 3 




V = 








/ 



Bi — A, 



Ji — \/ Ai-i/li 



From general theorems on the spectral theory of Jacobi matrices [T], there exists 
a positive measure d/i(x) such that the polynomials Xi( x ) are orthonormal with 
respect to it: 



(2.14) 



Xi(x)x 3 (x)d^(x) = 5. K 



Without extra conditions on the coefficients Ai and /li one could have more than 
one such measure, and correspondingly more than one random walk with different 
boundary conditions at infinity. Mathematically this boils down to the possible 
existence of several self-adjoint extensions of a given symmetric operator. Physically 
this corresponds to the possibility that a particle escapes to infinity in finite time 
in which case one needs to decide what to do at infinity. The reader can consult p] 
or [2] for a discussion of this issue. 



^The support of dfi(x) can in fact be shown to lie within [0, oo). 
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From (|2~14| it follows that 
(2.15) / Q l (x)Q J (x)d^(x) = -5 tJ . 



An integral representation for the transition probability functions is now obtained 
as follows. Introduce the functions 

OO 

(2.16) h{x,t) =Y,P ii {t)Q i {x) = [P(t)Q(x)]i. 
We have 

(2.17) = ^Q(x) = P(t)A(x)Q(x) = -xf(x, t). 
Given that 

(2.18) f(x,0) = Q(x), 

(2.19) fi(x,t) = e- xt Qi(x) 

and from the orthogonality relation (|2.15[) it follows that 

/OO 
e- xt Q l (x)Q 3 (x)d f i(x). 
-OO 

This is the celebrated Karlin-McGregor integral representation of the transition 
probabilities. In terms of the polynomials Xi{ x )i it reads 

(2.21) Pij ^)=i-ij J j- xt Xi (x) Xj (x)dn(x). 

At this point it is useful to think of a special case. If we had used the set of 
all integers as our state space, and we had chosen fa = Aj = 1 , so that the 
matrix above is the discrete Laplacian, the expression for Pij{t) would be given 
by exp(— 2t)li-j(2t), where Ifc(i) stands for the modified Bessel function. Since 
the eigenfunctions of the discrete Laplacian are given by the exponentials \k (0) = 
exp(ik8) the analog of (2.21) is nothing but the classical expression for Ifc(i) in the 
form 

(2.22) V( 27r ) / exp{2t cos 6) exp(ik9)d9. 

Jo 

In the quantum case , see for instance expression (13) in |12j one finds an ex- 
pression for the transition amplitude between sites m, / in terms of the usual Bessel 
function J7fc. Keeping in mind once again, the difference between the non- negative 
integers and the set of all integers as a state space, the classical expression for the 
Bessel functions is an instance of the expression p. 131) below. 

In the next section we shall examine how transition amplitudes for certain quan- 
tum spin chains compare with this. 



2 Here and in the following we trust the reader will not be confused by the occurence of i as 
label (index) or as the unit imaginary number 
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3. Quantum walks from spin chain Hamiltonians 

The continuous time quantum walk corresponding to the birth and death process 
characterized by the rate matrix A given in (|2.5[) is obtained by assigning vectors 
\i) to the states i and by using the matrix J of (|2. 12[) as the CTQW Hamiltonian. 
The action of J on the states \i) is naturally taken to be 

(3.1) J\i) = J i+1 \i + 1) + Bi\i) + Ji\i - 1). 

A central question is of course to determine the time-dependent transition ampli- 
tudesf ij (t) = (i\e iJt \j). 

We shall briefly show that this quantum walk can be identified with the single 
excitation dynamics of a spin chain with inhomogeneous couplings. This will make 
the connection between BDPs and spin chains. Secondly, we shall observe that the 
amplitudes /y (t) have an expression which is very similar to that of the probabilities 
Pij(t) given in (f2~2"T|) . 

Consider spin chain dynamics governed by Hamiltonians H of the XX type with 
nearest-neighbor interactions: 

N-l N 

(3-2) H=- ■fi+i(<W + i+°M + i)+2 EW + 1), 

1=0 1=0 

where Ji are the constants coupling the sites I — 1 and / and Bi are the strengths 
of the magnetic field of the sites I (I = 0, 1, . . . ). While the number of sites can be 
finite, we shall here allow for the chain to be possibly semi-infinite. One can adapt 
the formalism discussed here to the case of a doubly infinite chain following the 
ideas in [14] ,[19]. The operator H acts on C 2 ®°° (or C 2 ® N in the finite case) with 
the symbols erf, of , erf standing for the Pauli matrices which act as follows on the 
canonical basis {|0), |1)} for C 2 : 

(3.3a) <t x |1H|0), <7»|1> = i|0>, tr*|l) = |l), 

(3.3b) cr-jO) = ]X>, o*|0) = -t|l), ^ 2 |0) = -|0>. 

The index on these symbols indicate on which factor of C 2 they act. 
It is immediate to check that 

1 N 

(3-4) [H,- 5>f + 1)]=0 

1=0 

and it follows that the eigenstates of H split in subspaces labeled by the number 
of spins over the chain that are in state |1). 

Let us now restrict H to the eigensubspace spanned by the states with only one 
excitation, that is, states where all spins but one are in the |0) state . A natural 
basis for that subspace is given by the vectors 

(3.5) |i) = (0,0,...,l,0,0,...) ! n = 0,l,2,..., 

where the only " 1" occupies the i-th position. It is straightforward to see that in 
that basis, the restriction J of H to the one-excitation sector is given by the Jacobi 
matrix (|2.12l) . Moreover the action on the basis vectors \i) coincide with (|3.1I) . 
Note that the condition Jo = is assumed (and amounts to A_i = in (|2.13bp ). 
This hence establishes the desired identification with the quantum walk stemming 
from the BDP: the quantum states are the single excitation states of the spin chain 
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and the CTQW Hamiltonian J (related to the BDP rate matrix A) is simply the 
restriction of the spin chain Hamiltonian H to these states. 

The transition amplitudes are readily seen to satisfy a differential equation of 
Schrodinger type 



(3.6) 



ijJ i At)=i^\e- iJT \3) = ^\e- iJT J\3) 



Using a decomposition of the identity, one has 
(3.7) 



-J t hit) = E^le-^l^l J\j) = XXt)J« 



k=0 



k=0 



or in matrix form 
(3.8) 



df fJ 
"dt = fJ - 

Equations (|3 . T[) or (|3.8[) can now be solved exactly as was done for (|2 .3[) 
Let 



(3.9) 



Clearly, 
(3.10) 
We have 

(3.11) 



. d$j(x, t) 
dt 



$i(x,0) = Xi(x) 
.df(t) 



dt 



= [fJx]i = x$i(x,t). 



Hence 
(3.12) 

and using (12.14)) we find 



(3.13) 



/«(*) 



d/J-(x)xi(x)xj(x)e- 



We thus observe that, apart from a normalization factor, the transition amplitude 
fij (t) has an expression analogous to the Karlin-McGregor formula for the transition 
probabilities Pij(t) if we make the substitution t — > it. (This of course was to be 
expected given the definition of the CTQW). 

As observed after (|2.22[) an instance of this expression (|3.13|) is given in expression 
(13) in [12]. 

Note also that formulas similar to (|3.13p were obtained earlier in [S] and [21] in 
the context of Schrodinger quantum dynamics for the lattice Hamiltonians. 



4. Classes of spin chains with perfect and almost perfect return 

We shall illustrate the usefulness of formula (|3 . 1 3[) by considering the problem of 
perfect and almost perfect return for CTQWs generated by spin chain Hamiltonians 
of the form (|3.2[) . The probability amplitude for the walker to return to its initial 
state after some time T is fu(T). Such a return will be called perfect if \fu(T)\ = 1. 
An interesting question is to determine under which conditions can perfect return 
occur. 
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Let us first assume that at t = 0, the chain is in the state |0) = (1,0,0,...) and 
consider foo(t). From (|3.13|) we have 



oo 

ixt , 



(4.1) / 00 (t) = / e- lxt dfi(x) = F(-t), 

J — OO 

where F(t) is the characteristic function of the distribution defined as the 

Fourier transform of the measure [30] 

/oo 
e lxt d^{x). 
-oo 

The theory of characteristic functions is well developed [3D]. It is known in partic- 
ular that they possess the properties: (i) F(0) = 1 and (ii) \F(t)\ < 1, t ^ 0. In 
general \F(t)\ < 1 for all t ^ 0. However, there exists a special class of distributions 
dfi(x) for which \F(to)\ — 1 for some to > 0. This class only comprises (see [30] ) 
lattice distributions 

oo 

(4.3) fi(x)= M s S(x-t s ), 

S— — OQ 

where M s are arbitrary non-negative parameters such that Yl^L-oo ^* ~ 1 an< ^ 
where 

2ns 

(4.4) rs= £ + s = 0,±l,±2,... 

io 

with £ an arbitrary real parameter. 

Recall that the CTQW Hamiltonian is J, the one-excitation restriction of the 
XX spin chain Hamiltonian. We therefore have the following result 

Proposition 1. For a quantum walk initiated in the state |0), there is probability 
1 to return to that state after some time to > if and only if the measure d[i(x) 
corresponding to the Jacobi matrix J is a lattice measure |^73p. 

A similar result holds if the initial state is \i) with i = 1,2,..., i.e. if the initial 
state has its one spin up at the site i. In this case the returning amplitude is 



(4.5) A u = / xt(x)e- lxt dii{x) = / e' lxt d^{x), 

J —oo J —oo 

where 

(4.6) dfii(x) = x!(x)dfi(x). 

It is clear that d/j,i(x) is again a probability measure, because xf( x ) ^ an< i 

/>OG 

xl( x ) d K x ) = / dfii(x) = 1. 



in view of (|2.14l) . The amplitude Au(—t) is hence the characteristic function of the 
distribution dfii(x) and the requirement that dfi(x) be of the form (|4.3[) ensures also 
that a perfect return to any state \i) (\Aii\ = 1) is achieved. Note that modified 
measures of type (|4.6p (as well as the corresponding orthogonal polynomials) were 
studied in [18] . 

Consider now the generalization to the situation of almost perfect return. This 
means the following. Assume that the absolute value of the amplitude foo(t) is 
strictly less then 1 for all £ > 0. Suppose Nnevertheless that the value 1 can be 
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achieved with any prescribed accuracy for a sufficiently long time t. More exactly, 
this means that there exists a sequence of times t n , n = 0, 1, 2, . . . such that 

(4.7) lim |/oo(*»)| = l 

t n — >C30 

This is the condition for almost perfect return. 

Again, by a general result from the theory of characteristic functions [30] it 
follows that condition (|4.7[) is equivalent to the condition that the measure dfx(x) 
corresponds to a pure point spectrum, i.e. it has the expression (|4.3[) but now no 
other relations like (|4.4|) are assumed. 

Note that the amplitude foa{t) in this case is an almost periodic function |30j 
because it can be represented as a formal Fourier series 



(4.8) / 00 (t)= Ms 



which coincides with one of the possible definitions of almost periodic functions [4] . 
We thus have 

Proposition 2. A quantum walk initiated in the state 0} will almost perfectly 
return to that state if and only if the measure d[i(x) corresponding to the Jacobi 
matrix J has a pure point spectrum. This is equivalent to the statement that the 
amplitude foo(t) is an almost periodic function. 

5. Examples 

Let us first offer an example of CTQW corresponding to a linear BDP. Take for 
the birth and death rates the following expressions 

5.1 Xi = — ■ , m = - , 

1 — c 1 — c 

where j3 > and < c < 1. The corresponding symmetric Jacobi matrix J (|2.12[) 
has the entries 



(5 2) j. = Vct(i + (3~l) = (c+l)i + f3c 

1 — c 1 — c 

The Karlin-McGregor polynomials Qi (x) are readily identified with the Meixner 
polynomials Mi(x; ft, c) which are orthogonal with respect to the negative binomial 
distribution [29] 

(5.3) M(a;)= g ( i_ c) /3Mi£! 5(:c _ s) 

(as usual, the symbol (a) s = a(a + 1) . . . (a + s — 1) stands for the shifted factorial 
[29]). 

This distribution having discrete support, we expect the associated CTQW to 
exhibit perfect return. In this case, the characteristic function is obtained straight- 
forwardly 



(5.4) F(t) = I" e-^d^x) = (1 - cf p e- ist ( ^ (-c)* = ^ 

(-iy 



where 



10 



F. ALBERTO GRUNBAUM, LUC VINET, AND ALEXEI ZHEDANOV 



is the binomial coefficient. Expression (|5.4j) for the characteristic function of the 
negative binomial distribution is well known (see, e.g. (30]). Hence 

(5.5) / 00 (T) = F(-T) = 1 

ioTT = ±2nn, n = 0,1,2,.... 

Another example with perfect return corresponds to the so-called Stieltjes-Carlitz 
orthogonal polynomials related to the Jacobi elliptic functions (TU] . 

Consider two symmetric (i.e. with B n — 0) orthogonal polynomials C n (x) and 
D n (x) defined by the recurrence relations 

(5.6) Jn+\Pn+\{x) + J n P n -i(x) = xP n (x), 

where P n (x) stands for either C n (x) or D n (x) and the recurrence coefficients are 
defined as 



(5.7) J 2 n = 
for the polynomials C n (x) and 

(5.8) Jl = 



z n if n even 
n 2 if n odd 



k 2 n 2 if n odd 



for the polynomials D n (x). In these formulas k is an arbitrary real parameter 
satisfying the condition < k < 1. 

It can be shown that the polynomials C n (x) are orthogonal on the infinite discrete 
set of points 

(5.9) Ts = ^(2s + l)m s = 0,±l,±2,..., 
where 

(5.10) K ' ^ 



\J\-k 2 sin 2 <p 

is the complete elliptic integral of the first kind. 

These polynomials satisfy the orthogonality relation 



oo 



(5.11) ^ M s C n {T s )C m {T s ) = S nm 

S — — QO 

with 

where q = ex~p(—TrK'/K) is a standard parameter in the theory of Jacobi elliptic 
functions and r\c is an appropriate normalization constant. 

Similarly, the polynomials D n (x) are orthogonal on the infinite discrete set of 
points 

(5.13) Ts =x' s = 0,±l,±2,... 
and satisfy the orthogonality relation 

oo 

(5.14) M s D n {T s )D m {r s ) = S nm 
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with 

(5.15) M s = —^7- 

q s + q * 

In both cases the spectrum is discrete and linear with respect to s. Hence both 
polynomials C n (x) and D n (x) entail condition of the perfect return. Note however 
that in contrast to the Karlin-McGregor theory, the orthogonality interval is the 
whole real axis [— oo, oo]. This means that the recurrence coefficients of the Sticltjcs- 
Carlitz polynomials C n (x), D n {x) do not meet the conditions needed in order to 
define a Markov birth and death process. For the quantum walks, however, such a 
restriction on the orthogonality domain is not needed and we have here a perfectly 
valid example with a double-infinite spectrum. This difference between the classical 
and quantum walks is worth noting. 
It is easily seen that the expressions 

oo 

(5.16) /oo(f)= M s e~ itTa 

S — — OC 

coincide with the Fourier series of the elliptic Jacobi functions cn(z; k) and dn(z; k) 

m 

(5.17) /oo(i)=cn(wi;fc) 
for the polynomials C n (x) and 

(5.18) /oo(t)=dn(wt;fc) 

for the polynomials D n (x) with some parameter ui. We see that indeed, both func- 
tions are periodic and achieve values 1, as should be for the perfect return. There 
is an interesting difference between these two amplitudes. Indeed, the amplitude 
(15.171) can take zero values for some values of the time t, while the amplitude (I5.18[) 
never achieves zero values on the real axis t. This means that the amplitude (j5 . 1 7[) 
possesses the property of "self-avoiding" (i.e. foo(U) = on a uniform lattice ij). 

Spin chains with a finite number of sites will have a discrete spectrum and hence 
will always exhibit almost perfect return. It is also clear that such chains with 
perfect state transfer (PST) [3J, [33] will always manifest perfect return by time 
reversal symmetry. 

As we already know, for almost perfect return it is necessary and sufficient that 
the orthogonality measure dfj,(x) be discrete. In all other cases almost perfect return 
is impossible. In this connection, consider a semi-infinite uniform Heisenberg XX 
chain without external magnetic field; this means that Bi = and that interaction 
constants are all the same Ji = 1/2, i = 1,2,.... The orthogonal polynomials 
Xn{x) corresponding to the Jacobi matrix J coincide in this case with the Chebyshev 
polynomials of the second kind: 

(5.19) x»(*) = Slng( " + 1) , * = cos0. 

suit* 

It is well known [TU] that these Chebyshev polynomials are orthogonal with respect 
to the purely continuous measure 

ri 

(5.20) J Xn{x)Xrn{x)\/l - X 2 dx = —6 nm . 

We thus see that for the uniform Heisenberg chain even almost perfect transfer is 
impossible. 
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This can be confirmed by direct calculation. Indeed, we have 



(5.21) /„,,(*) = - / - ^e~ ixt dx = 



* J-i 




where J\(t) is the first Bcsscl function of the first kind. The function foo(t) is real, 
takes value 1 at t = and oscillates with a decreasing amplitude for t > 0. The 
first minimum of foo(t) occurs for t\ rts 5.14. Its absolute value is |/oo(*i)| ~ 0.13. 
All further local extrema £2,^3, ■ • • correspond to smaller values of |/oo(i)|- This 
means that the homogeneous Heisenberg chain is far from showing almost perfect 
return. 



We trust this communication is providing an interesting bridge between birth and 
death processes and quantum walks generated by XX spin chains. The relation 
between CTRWs on graph and CTQWs is well established in general. We wished 
to draw the parallel further in the special case where the graph is linear and the 
generator matrix tri-diagonal, that is when we have a BDP. In this instance, the 
CTQW can be viewed as the single-excitation dynamics of spin chains with nearest- 
neighbor interaction. Furthermore, as shown by Karlin and McGregor in their 
analysis of BDPs, the theory of orthogonal polynomials offers powerful tools when 
Jacobi matrices are at play. We here have given indications on how this approach 
carries in the quantum framework. The study of BDPs is a classical topic, we hope 
that the correspondence with CTQWs will prove helpful in the transfer of results to 
the quantum walk realm. The spectral results have been put to use to characterize 
the systems for which the probability to return to the starting point is 1 or almost 
1. Almost perfect return occurs if the 1-excitations have a pure point spectrum 
while perfectness is achieved if the spectrum is uniform in addition. 
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